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^ We study the singularities in X-ray absorption spectra of one- dimensional Hubbard and t-J 

models. We use Boundary Conformal Field Theory and the Bethe Ansatz solutions of these models 
with both periodic and open boundary conditions to calculate the exponents describing the power- 

VQ ! Is-w decay near the edges of X-ray absorption spectra in the case where the core-hole potential has 

rvq . bound states. 

^ I. INTRODUCTION 



B 

I 



X-ray absorption in a metal can be described by a simple model put forward by Nozieres and de Dominicistl. An 
electron from a filled inner shell of one of the nuclei is raised into the conduction band. This generates a local potential 

'"^ , V at the position of the nucleus that lost the core-electron, which in turn acts on the (noninteracting) conduction-band 

^ ' electrons and affects the X-ray absorption probability. The situation is described by the Hamiltonian 

O' 

,'^, . i7 = ^e(fc) ct(fc)c(fc) + 66t^y(fc,fc') c^(^)c(^')+£^o&^^>, (1-1) 

(N : 

*^ ] where e(fc) is the dispersion of the conduction band electrons, h^ and b {c{k) and ct(fc)) are annihilation and creation 

r-O ' operators for the core-hole (for conduction band electrons with wavevector k) and Eq is the energy of the core state. 

psl ' As b^b commutes with H, the Hilbert space splits into two sectors: in one the core-level is filled {bb^ — 0) and there is 

CN no potential whereas in the other one the core level is empty (66^ = i) and V acts on the conduction electrons. As was 

^^ ' shown in Ref. ^ the inner core disturbance acts as a transient one-body potential on the conduction electrons, which 

^^ , means that one needs to study the response of the conduction band electrons to the potential V applied between 

^~/ ■ times t = and t ~ t' . The X-ray absorption rate can be expressed by the golden rule as 
-)— > I 

^ I{u)^y^\{n\clmQ)? S{u + Egs-E^-Eo), (1.2) 



a 

(-H where co(i) annihilates a conduction-band electron at position a? = at time t, |0) is the ground state at times 

Q t < and H\Q) — Egs\^)- The r.h.s. of ( pT^ ) can be expressed in terms of the spectral representation of the Fourier 

O transform of the retarded correlation function ((&''' (i)co(i)cJ (0)6(0))) so that 



^ 



'k^ ; I{uj) ex 5m / dt e*"*((6^(i)co(i)cJ(0)6(0))) . (1.3) 



Near the threshold wq ~ E^ the intensity I{lo) displays a characteristic singularity of the form 

/(a>) ~ I ^-— . (1.4) 



For the system ( 1 . 1 )|-jl|iie critical exponent a has been determined exactly and is expressed in terms of the phase shift 
at the Fermi surfacma. A very interesting case is the one where the local potential V is sufficiently strong to bind a 
conduction electronEl (see also Refs. ^J^)- ^^ ^^^^ '^^^^ ^^^ absorption spectrum (if a > 0) features two thresholds with 
characteristic power-law decays of / as a function of uj (see Fig. ^). If a < there is no discontinuity and I{uj) goes 
to zero instead (see Fig. mb). 
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FIG. 1. X-ray absorption rate as a function of frequency: a) qq > 0, Qi > , b) ao < 0, ai > 0, 



In the present work we wish to investigate the analogous situation for integrable lattice models of strongly interacting 
conduction electrons in one dimensioiJj'Q. These models are particular realizations of Luttinger liquids and the X-ray 
problem for such systems has been investigated by various authors (a detailed pedagocial discussion can be found in 
the forthcoming boom). The case of a core potential with no backscattering was solved in Refs. g and the case pi 
a perfectly reflecting potential was treated in Ref_l^. The general case was investigated by Affleck and Ludwigtll 
using Bou ndary Conformal Field Theory (BCFT)t37llecently, Afflecka reconsidered the X-ray problem for a Fermi 
liquid (1.1) for the case where V has a bound state from the point of view of BCFT. This motivated the present work 
in which we study the X-ray problem in Hubbard and t-J chains for core hole potentials with bound states. Let us 
discuss the general setup for the case of the Hubbard model. At times < < we take the system to be periodic 



Ha = 



EEH 



cr'^j+l,a ~r C^^-^^^Cj^^ 



iu 



J = l 



L 



liN. 



(1.5) 



At time i = we switch on the core potential Vil acting on sites 1 and L (a similar situation has been studied inll2 
In the general cast; this potential will include a backscattering term which will then drive the system to the open 
chain fixed pointlla, i.e. break the chain across the link IL. Wc model this situation by considering the Hamiltonian 



L-l 



Hb = - 



E E (4.-^JH 



'-J + l.cr'-J.O- 



4-E 



rijpiji + ^N + Hi + Hl 



(1.6) 



where ffi,L are one-body interactions acting on sites 1 and L respectively. At time t we switch off the core potential 
which changes the Hamiltonian back to Ha- Depending on the precise form of the interactions Hil bound states can 
be formed at the boundaries. As the elementary excitations in the Hubbard model are not electrons like in the case 
of the Fermi liquid discussed above but (anti)holons and spinous one has to consider several possibilities: In addition 
to the case in which there are no bound states the core-hole potential can bind either a spinon, a (anti)holon, both a 
spinon and a (anti)holon or, for an attractive boundary potential of the order of the Hubbard interaction 4u, a pair 
of electrons. 

In order to extract the X-ray exponent we use BCFT and the fact that the low-energy spectrum of both Hubbard and 
t-J models can he. described in terms of two c = 1 Conformal Field Theories or equivalently a spin and charge separated 
Luttinger liquidEflO. Our discussion closely follows Ref. |ll|. We start by considering the Luttinger liquid defined 
on the complex plane with coordinate z. Identifying the radial part of z with the time variable the case of periodic 
boundary conditions (A) is realized if we consider the complex plane without boundaries. The change to open boundary 
conditions (B) corresponds to the introduction of a cut in the plane from zq to zi. As explained above this change of 
boundary conditions corresponds to switching on (and off) the core-hole potential. Choosing < tq = zq < ti = zi 



real and mapping the plane to a cylinder via the conformal transformation z = exp(27r(u + iv)/L) this cut gets 
mapped onto a seam in the time direction of the cylinder (see Fig. y). 
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FIG. 2. The mapping from the cylinder to the plane. 
The Green's function of an operator O with dimension x on the complex plane without boundaries is given by 



{A\0{T,)0Hr2)\A) = 



1 



(ri - T2)2- ■ 
The Green's function on the cylinder is obtained by the conformal mapping. For U2 ^ ui ':$> L we obtain 



{A\0{u,)0^{u2)\A) 



-2^(m2-ui 



(1.7) 



(1.8) 



To study the edge singularity we choose O^ to be an operator which changes the boundary conditions from A to B. 
The same correlation function can be evaluated alternatively by inserting a resolution of the identity in terms of the 
eigenstates \B;i/) of the system with reflecting boundary conditions 



{A\0{m)OHu2)\A)=J2momB;,.)\^ e"'^--^- 



)(m2-«i) 



(1.9) 



The leading contribution to this sum comes from the ground state or a low lying excited state (this depends on the 
operator O because the form factor must be nonvanishing) with boundary condition of type B. Comparing the two 
expressions for the correlation functions on the cylinder allows one to extract the scaling dimensions of the boundary 
changing operator O 



L 

2^ 



{Eg - E^) , 



(1.10) 



For boundary potentials that do not lead to bound states one identifies the exponents Xq for the core-hole operator 
and xi for the core-hole conduction-ielectron operator {E^ be ing the ground state energi es in the -/V-((-/V-|- l)-)particle 
sector with B boundary conditions )llil. Fourier transforming (1.7) the edge exponent in (|l.4D is identified as 



a = 1 — 2 a;i 



(1.11) 



In the presence of the various types of bound states the power-law behaviour (|l.4| ) oi I{lu) above the respective 
thresholds can be determined by inserting the appropriate excited-state energy into ( 1.10| ). Finally, let us note that 



in the above discussion we have set the Fermi velocities to one; the generahzation to the two-component Luttinger 
Uquid with different Fermi velocities proceeds along the same lines as in the case of periodic boundary conditionsEZIEj. 
In the remainder of the paper we follow the steps outlined above to study the nature of the X-ray edge singularities 
in the t-J and Hubbard models for boundary terms Hil chosen in such a way that they preserve the integrability of 
these systems. 

II. THE T-J MODEL 

In this section we determine the X-ray absorption exponents for|*|i-J chain with the particular choice of core- hole 
potential described above. We consider the following Hamiltonianala 

L-l L-1 

+2 J2 Sj ■ Sj+i - "'"'+' + 5] n, + fij+i -fiN + H^p , (2.1) 

j=i j=i 

where V projects out double occupancies, 5*^ are spin operators at site j, rij = c- tCj^^+c, yCj^i and Uj = 1 — n^j — nj_|. 
There are three different forms for the boundary part Ha/3 of the Hamiltonian that are compatible with integrability 
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H,, = him + hLUL , H^s. = him + hL{Sl - ^) , H^^^h{Sl-^ + Sl-^) . (2.2) 

These correspond to localized potential (a) and magnetic (b) interactions of the conduction electrons with the distur- 
bance due to the core hole. Physically local magnetic field interactions are not very realistic; one would rather expect 
a Kondo-like interaction which we cannot consider in the present framework of integrable lattice models. In what 
follows we therefore constrain our analysis to the model with aa boundary conditions. We note that in the continuum 
limit iJaa gives rise to forward scattering terms. We therefore expect that the X-ray exponents will generally not be 
universal despite the fact that our boundary is perfectly reflecting in the sense of Refs. |o[^ As we will se below 
this JiLmdeed the case. However, the situation is somewhat more complicated than this: unlike in the continuum 
limitlljO we do not impose Neumann boundary conditions (on the lattice wave functions) . The boundary conditions 
should rather be thought of as being of mixed Dirichlet-Neumann type (e.g. c7/;(0) -I- dxtpi^) = 0). The parameter c 
enters the finite-size spectrum in the same way as the forward scattering amplitude. Therefore in the continuum limit 
the forward scattering amplitude is not simply given by the boundary chemical potential. As a result we recover the 
results of Refs. pX)| , pT| not for hi^L ^^ 0, 1 but for some finite value that depends on the band filling (see below). 

In the following we start by considering boundary fields in the region 1 < /i < 2. This is unphysical from the 
point of view of the X-ray edge singularity where the potential due to the core hole should be attractive but permits 
a pedagogical discussion of the formalism we use to calculate the finite-size energies necessary for extracting X-ray 
exponents. 

A. Repulsive boundary fields: 1 < /i < 2 

In this region of boundary fields holon boundary bound states at both boundaries are present in the ground state 
of the t-J chain. Defining 

5-, =2- A, j = l,L (2.3) 

the Bethe Ansatz equations with respect to the reference state with all spins up readli3 

(ei(A„))'^ = \{ e2(A„-A^)e2(A„ + A^)[|e_i(A„-AW)e„i(A„ + A«) 

/3/a 7=1 

l = e_5.(A«)e_s.(AW) [] ^^(AW - A^)ei(AW + A^) , (2.4) 

13=1 



where e„(x) = ^_^ . The energy of a state corresponding to a solution of (2.4) is 



Nh + Ni 

E^hi + hL- V -^-—^ + fiNh . (2.5) 



U I + A| 



We now observe that for hi > 1 solutions of (2.4) exists where (in the thermodynamic limit) two roots A^^-' take the 
values —^Si and —^Sl respectively. These roots correspond to boundary bound states. The situation is analogous 
to the XXZ Heisenberg chain studied in Rcf. M. One finds that the ground state is given by a distribution of roots 



such that both these boundary roots are present. The logarithmic form of the Bethe equations (for a solution of (2.4) 
with only real roots apart from the boundary roots) reads 



^/^ = (2 + y)^(A.)-iE^(^^) + ' 



,Xa + A/3 , 



L' ' ' L^-^ ' 2 ' ' 2 



7=1 

T^7 = zE^(^7'^-^")+^(^7'^+^") + ^^^' 7=1.. .^/.-l (2.6) 

a 

where N^ is the number of electrons with spin down, Nh is the number of holes, /^'"^ are integer numbers, 0{x) = 
2 arctan(2x) and 

"(0 = -«(^) - D(^) ■ (2-7) 



In addition to (2.6) we still have two equations determining the precise values of the boundary roots. A detailed analysis 
of these equations yields that the corrections to the thermodynamic values in a finite system vanish exponentially 
with system size. This mea ns t hat for the purposes of the present work we can neglect these corrections. We should 



note here that solutions of (2^) do not yield a complete set of states_For vanishing boundary fields such a basis can 



be constructed by means of the s/(l|2) symmetry of the HamiltonianES. For nonzero boundary field s this symmetry is 



broken and we do not known how to complement the set of Bethe states given by solutions of (2.6). However, for the 
present purposes this is not necessary: we are interested in the lowest energy state in a p articular sector of quantum 
numbers and it can be shown that these states can always be obtained as solutions of ( |2.6| ) or the analogous equations 
based on the Bethe Ansatz reference state with all spins down. We note that this ceases to be true for the t-J chain 
with ba or bb boundary terms. 

The calculation of the finite-size spectrum proceeds along the lines of Refs. Il^ and pO so that we merely quote the 
result 



-'•"— /-^{^^^^V^-i-^} 



+^l^fAN;-^-ei + &f-± + Nll (2-8) 

where Coc is the ground state energy of the infinite system, Vc and Vg are the Fermi velocities of holons and spinous 
respectively, ^ = ^(A^) is the dressed charge defined via 

e(A) = l+/ d,,Gi{X-,y)^{iy) , 

^ .^^ 1 1"°° J -^c.A e-l^^l 1^ f,/3 + 2; .A, ,,1 + a; .A.I 

GAX) = ^J_J^e -^^^;^^ = ^5ie|^(^ + z-)-V^(^+^-)j, (2.9) 

where ip{x) is the digamma function. The integration boundary Ac is determined by the chemical potential (band 
filling). We note that as we approach half-filling (/i -^ 2 In 2) Ac ~ \/ 3^3) (2 In 2 — /i). 



The term proportional to N" = ^^^ jj-g Ip — /^ is the contribution of particle-hole excitations, where /"^ are the 
integers corresponding to the roots of the particle and the hole. The quantities A7V° and AA''^ denote the deviations 
of the total particle number and the number of down spins from their respective values for some reference state. This 
concept needs to be introduced because in order to extract the X-ray exponents we need to compare finite-size energies 
for different boundary conditions. The state with respect to which we measure the deviations of particle numbers is 
chosen such that for the ground state A7V° — Oq — for a = c,s. This may appear odd but turns out to be the most 
convenient choice for calculating the energy difference between states with open and closed boundary conditions. The 
quantities 9^^ are defined as 

pi(A)=.g„.o(A)+ / dv[5^sG^{\-v)+5^,Gi{\-v)]pl{v) , a = c,s, (2.10) 

J-Aa 

where 

<7.,o(A)= J2 Gs,(A)+G|i_5,|-i+Gi(A), 

ffe,o(A)= Y. Gi+5,(A)+G|i_5,|-«s,(A)-Go(A) . (2.11) 

Last but not least the surface energy foa is found to be 

foo=fo + fcihl) + fc{hL) = fo + /bound , (2-12) 

where /o is the surface energy of the system in the absence of the boundary bound stateaia and fc{hj) are the 
contributions of the holon boundary bound states. Note that these contributions are of order one unless we fine-tune 
the boundary fields. We findlia 

1 /■^'= 1 

/o = -2 / dAec(A)[asi(A)+a5jA)]--[£,(0) + /i-2/ii-2;iL] , 

UK) = /i - ^(G3_ 2 (0) + G_i+ . (0)) + i j^^^ dv [G3_ 2 (,.) + G_i+ . {v)] e,{v) , (2.13) 

where ax{X) — -^ ^ ^^^ and where the dressed energies are given as solutions of 

rA. 

es{X) = -2^Go(A) + / di^ Go(A - j^) e,{i^) , 

£c(A) =/i-27rGi(A)+ / dv Gi{X - ly) ec{v) ■ (2.14) 

J-Ae 

The bound state energy fc{h) as a function of boundary chemical potential is shown for different band fillings in 
Fig. g. This characterizes the relevant part of the low-lying finite-size spectrum of the open t-J chain with boundary 
fields in the sector where N^^ > N^. In order to extract the X-ray exponents we need to consider states with N^ > N^ 
for the case where the core electron carries spin down. This can be taken care of by changing the reference state of 
the Bethe Ansatz to the state with all spins downta. The result is of the same form as ( p.8[ ) but with redefined AN^. 
We also need the finite-size ground state energy of the t-J model with periodic boundary conditions. It is given 

byta 

E^°^ =Leoo-^ivc + v,) . (2.15) 

bL 

We now have the necessary machinery to determine X-ray exponents. One should keep in mind that we presently 
have repulsive boundary fields. For pedagogical reason we nonetheless will calculate X-ray exponents for this case: 

Absolute Threshold: The lowest (in frequency) threshold in the X-ray absorption intensity occurs at some 
frequency ujq and is associated with an intermediate state in which both holon bound states are occupied. For the 
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FIG. 3. Energy of holon bound states as a function of boundary chemical potential. 



case wher e the core electron has spin up this corresponds to AN^ ~ —3 , AA^" = —1. Combining ( ^.8[ ), (2.15), (I.IC) 
and ( 1.11 ) we obtain 



O^abs ^ 



I [3 + 08] 



,cl2 



2^2 



(2.16) 



For the case where the core electron has spin down we need to proceed as outlin ed above and use the Bethe Ansatz 
solution with a different reference state. The final result is the same as ( 2.16| ) as TJaa preserves the discrete spin 
reversal symmetry. 

Intermediate Thresholds: The second and third thresholds occur when one of the holon bound states is 
occupied but the other one is not. Let us consider the case where the bound state at 1 is occupied. The corresponding 
threshold in the X-ray absorption rate is at luq — fc{hi). As only one holon bound state is occupied we now have 
A7V° = —1, AN^ = —2 and the expressions for the quantities 9'^''' in ( |2.§| ) get modified. They are now given by ( 2.10| ) 
but with different driving terms 



gc.iW - Gi+5i(A) + Gi_5i(A) - Go(A) - as. 



The X-ray exponent associated with this threshold is 



^int 



[2 + 01 
2^2 



cl2 



(2.17) 



(2.18) 



Band Threshold: The fourth and final threshold occurs at luq — fc{hL) ~ fc{hi) when neither bound state is 
occupied. This corresponds to the case where the core electron is emitted into the conduction band where it decomposes 
into an antiholon and a spinon. Then AA^^^ = -1 , AiV^" = -1, e'§ = 26*1, g^sW = -Go(A) - as,{X) - as^(A) and 
the associated X-ray exponent is 



^band 



c\2 



(1 + 1) 

4^2 



(2.19) 



B. Attractive Boundary fields: < h < 1 



This r egio n of boundary fields corresponds to an attractive core hole potential because of the form of the third last 
term in (2.1). Now no boundary bound states exist. The analysis of the finite-size spectrum follows the one a bove , 
the only difference being the absence of purely imaginary roots. The X-ray exponent is of the same form as ( 2.19| ) 
where we should keep in mind that Si^l now are negative. The results for two different band fillings are plotted in 
Fig. ^ as functions of the boundary chemical potential h = hi = Hl- Our result coincides with Rcfs. [lO| , [ll| if we make 
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FIG. 4. X-ray exponents (a) and vr^*^/^ (b) for the t-J model with aa boundary conditions and hi — Hl = h < 1. 



the identification 6'^ = — ^^, where Vf is the forward scattering amplitude of the core hole potential in the continuum 
limit. We see that 9'^ does not vanish for hi^ —* 0. As explained above the continuum Vf is not simply given by the 
boundary chemical potential so that there is no contradiction. In Fig. Hb we plot tt^^/^^ as a function of h. 



C. Attractive Boundary fields: h < 



In this range of boundary chemical potential the analysis of the finite-size spectrum is less intuitive than above. 
The Bethe equations (2.4) allow a variety of boundary string solutions like the ones encountered in the repulsive case. 
However one finds that none of these complex roots is present in the ground state. We interpret this as follows: in 
the ground state antiholons and spinous are bound to the boundaries. States where some of these bound states are 
unoccupied are characterized by imaginary roots of the Bethe equations. In support of this interpretation we can 
compute the particle number at the boundary site. It is given by ^, where h is the boundary field. We find that in 
the ground state there is a strong enhancement of charge at the boundary site as compared to the bulk. The states 
involving imaginary roots of the Bethe equations exhibit a significant decrease in charge at the boundary as compared 
to the ground state, which is consistent with our interpretation. 

Absolute Threshold: 



In order to calculate the X-ray exponent for the lowest threshold we need the finite-size energy of the ground state 
for h < 0. As no complex roots of the Bethe equations are present the analysis is straightforward and very similar to 
the band threshold for 2 > ft, > 1. We find 



^abs i~. 



1 {l + 9'=f 



2e 



(2.20) 



where 9'^ is given by ( ^.10 ) with gc(A) = — Go(A) — asi(A) — aSi,(A). 

In Fig. the X-ray exponents of the absolute threshold are plotted as functions of the boundary chemical potential 
for two different band fillings. For simplicity we only consider the case hi = hj^ = h. We see that in the physical 
regime h < 1 there is always a singularity associated with the absolute threshold i.e. I{lu) always diverges. 

Higher Thresholds: Let us consider the case in which two complex roots A*-^^ are present and take the values 
— ^Si^L respectively. The Bethe equations read 
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FIG. 5. X-ray exponents for the absolute threshold in the t-J model with hi — Hl = h at almost half-filling and quarter-filling. 



{ei{Xa)f Yi ei_s,(AQ)ei+s.(Aa) 

= n e2(A„-A^)e2(A„+A^) n e_i(A„-AW)e_i(A„+AW) 

fJ^a 7=1 

l = e-s,(A«)e_5.(AW) H ei(AW „ A^)ei(AW + A^) . 

/3=1 



(2.21) 



Following through the same steps as before we find that this state has a gap of magnitude A/ = fc{Si) + fc{SL) 
where 



f,iS) = Ai - 7r[Gs+i(0) - Gs-iiO)] + ]- T d\ e,{\)[Gs+i{X) - Gs-i(A)] 



(2.22) 



We interpret this state as differing from the ground state by leaving boundary bound states of antiholons unoccupied. 
Consequently we find a threshold in the X-ray absorption probability at a frequency A/ higher than the absolute 
threshold with exponent 



*-^int (., 



1 (3-Hr) 



c^2 



2^2 



(2.23) 



where ^ is given by (^ with g,{\) = -Go(A) + Ej=i,l Gs,+i(A) - Gs,_i(A) - as, (A). 

Thresholds at lower frequencies occur if we have only one imaginary root A^^^ — —-^S where 5* is either Si or Sl- 
The corresponding states have a gap equal to A/ = fdS) and give rise to exponents 



a„ 



{2 + e^ 

2e 



(2.24) 



where 9'^ is given by ( |2.10| ) with gc(A) = —Gq{\) + Gs+i{x) — Gs~i(A) — as^(A) — as^(A). A numerical solution of the 
relevant integral equations for a quarter filled band shows that aint is negative and therefore leads to a "shoulder" in 
I{uj) as in Fig. |l| b). On the other hand we find that a'^^^ is positive and leads to a singularity. 



The cases imvestigated above by no means exhaust the hst of states with imaginary roots. For example there is 
a state with two imaginary A^^^'s taking the values —^Si^l and two imaginary A's taking the values |(1 — Si^l) 
respectively. This type of solution of the Bethe equation also gives rise to three thresholds as imaginary A's are only 
allowed if their respective "partner" A^^^ is present as well. The calculation of the X-ray exponents is completely 
analogous to the case treated above so that we omit it. 



III. THE HUBBARD MODEL 



The one dimensional Hubbard model with open boundary conditions of type aa (i.e. boundary chemical potentials 
only) 
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H = - 



T.T.H 



o-S + l,'^ 



S'+1,ct'^J'' 



4m 



j=i 



3 = 1 



T-jlTiji + iiN — hiTLi — h^riL 



(3.1) 



is soluble by means of the Bethe Ansatz as sho wn in Refs. PlJ,E3 (note that the boundary potentials are defined in a 
different way than above: to identify hi^L in (3.1) with those used for the t-J model one should replace hi^^ —^ 1— ^i,l)- 
Applying boundary magnetic fields instead also leaves the Hubb ard model integrablec^ but will not be considered 
here. The Bethe Ar^s^itz. equations determining the spectrum of (|3.l[) in the A'^e-particle sector with magnetization 
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a = 1, , 



,^i. 



(3.2) 



The quasi momenta kj and the spin rapidities Aq, paramatrize an eigenstate of (3.1) with energy 



E = nN, 



2 y^ cos kj 



(3.3) 



For small values of the boundary fields the ground state configuration is given by distributions of real kj and Aq and 



Bi^\k) 



1 - h.^e' 



i?(-)(A) = 1 , 



(3.4) 



contain the phase shifts due to the boundaries (this case has bee n di scussed in Ref. ^2|). For sufhciently large boundary 
chemical potentials /ii,l, however, the Bethe Ansatz equations (3.2) allow for various complex solutions corresponding 
to boundary bound states for antiholons, spinous and pairs of electrons, respectivclyEj: First, for hix > 1 one finds 
bound states parametrized by k = i In/ii^^ with exponential accuracy in the thermodynamic limit L — )■ oo. The quasi 
momenta parametrize the charge part of the states: hence this solution corresponds to a charge (or antiholon) bound 
to the surface. Inserting this solution in the second set of Eqs. (^) leads to a boundary phase shift in addition to 
the product over the real quasi momenta kj which modifies B^ {Be remains unchanged): 

Si"nA) = e2(«+5.)(A)e2(„-s.)(A), (3.5) 

where we have introduced Sx = {hx — l/hx)/^, > with .x — 1 or L. Analyzing the resulting equations we find that 
a new type of solution arises at Sx — u, i.e. hx — u + \/v? + 1: Beyond this point a complex solution A — i{Sx — u) 
for the spin rapidities is allowed. We note that spinous are to be identified with holes in the distribution of spin 
rapidities. Again, occupation of this state modifies the boundary phase shifts Bc,s- 
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Bi'^Hk) = L _ ^ Jk ) e-2sJsinA:j)e2(s,_2«)(sinfcj) , 

Bi^\X) = e2(s._3„)(A) e2(„_s.)(A) . (3.6) 

Finally, boundary potentials with Sx > 2u can bind a (singlet) pair of electrons to site x. Such a state is parametrized 

by two complex quasi momenta sinfcg = Xq zLiu and a single complex sp in r apidity Aq = i{Sx — u) as before. The 
remaining real solutions of the Bethe Ansatz equations are determined by (3.2) with 

^'^^^^'^ ^ ij-h^) e-2S.(sinA;,)e2(5_2n)(sinfc,) , bI^HA) = 1 • (3.7) 

Depending on the strength of the boundary potential we have to distinguish between the following cases in order 
to describe the spectrum: in addition to the case discussed in Ref. ^, where the solution of the Bethe Ansatz eqs. is 
given in terms of real kj and Aq only, one can find either 

• an antiholon in a bound state (corresponding to a complex k) and the spinon in the corresponding band (which 
implies the presence of a complex A for Sx > u), 

• an antiholon and a spinon in bound states (parametrized by a complex k for Sx > u), 

• and finally, for Sx > 2u, a pair of electrons bound by the potential. 

Each of these configurations gives rise to a continuous spectrum above a threshold that depend on the occupation of 
the boundary states. 

In the following, we shall discuss some of these cases for the symmetric choice hi = Hl = h oi the boundary 
potentials. The bound states discussed above will occur pairwise at the given thresholds (corresponding to sit es 1 and 
L, respectively). As for the t-J model we shall consider the logarithmic form of the Bethe Ansatz equations ( |3.2| ) for 
low lying states above these thresholds: 

27r/,- 1 v-^ f / sin fc, — Ag \ / sin fc, + Ag \ 1 1 , , 

^ = 2., + -g|.(^-^j+.(^-^j| + -.(/.,), J-l,-,N 

27: Ja 1 v-^ f „ / Aq — sin kj\ /, / A^ + sin kj 



E 



L L ^^ I V 2u I \ 2u 



1 ^ 



>:-^k'\ , .^^.. + A.A\ , 1 (^_,) _ „ = 1,...,M. (3.8) 



L '^^ [ \ 4:U J v4u/J L 

Here the summations extend over the real roots kj and Aq,. The functions k, and lo contain the phase shifts due to 
the boundary fields and occupation of the boundary bound states. 

A. Band threshold 

The edge singularity with the highest threshold corresponds to excitation in states with no bound states occupied 
by the particles. This situation was studied in Ref. p2l Like in the case of the t-J model this does in fact imply the 
occupation of a holon bound state for repulsive boundary potentials hij-^ —1: computation of the particle number 
on the boundary site shows a depletion due to the presence of the h olorO . In the Bethe Ansatz equ ations the only 



boundary phase shifts are those due to the boundary potentials, i.e. (3.4). The resulting K{k) in ( |3.8| ) is given by 



while Lu{X) — 0. The finite size spectrum for the relevant boundary conditions is again given by (2.8) and (2.15) 
The dressed charge i^ = £,{Q) for the Hubbard model is defined in terms of the solution ot-^^ integral equation {Q 
varies between and tt as a function of the density of electrons and the coupling constant )c 
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^(fc) = 1 + / dk' cos k' K{sm k - sin k')^{k' 
J-Q 



'-Q 

K(x) = — / duj 

2n Jq cosh uui 



Here 6''^'^ are related via 9^ = ^6'^ with 



\[ r dkp^Pik) 1 



(3.10) 



(3.11) 



-Q 



.(1) 



for our choice of the reference state. The 0{1/L) contribution pi to the density from the boundary fields is given in 
terms of the integral equation 



p)^> (k) = 5c(fc) + cos fc / dk'K{sm k - sin k')pi^> (fc') 



(3.12) 

For the case considered here the driving term in this equation is found to be (after integrating out the spinon-part of 
the densities) 



= (fc)=#(fc) 



1 



1 



cosk 



TT 1 + ft,2 _ 2/i cos fc 4-ucosh(2^ sinfc) 



(3.13) 



An analytic solution of this integral equation is possible in certain limits only. It simplifies essentially in the strong 
coupling limit where K{x) = ln2/27ru. This allows to give a simpler expression for 9'^ in terms of the driving term 



1 ln2 . ^ 
1 H sm Q 

TTU 



dkgc{k) — 1 



for 



Furthermore it is known that Q — nric and ^ = 1 in this limit. With (3.13) we find 

1 
2 



2 /I 

9^ = — arctan — 
TT VI 



h nn^ 

— tan 

h 2 



(3.14) 



(3.15) 



for infinite couplin 



In general the integral equations have to be solved numerically to compute the X-ray edge 
to the finite size ground state energy of the Hubbard chain with periodic 
For absorption of the core electron into the band we have to choose AA^^ — 1. The 



exponents from (1.10) by c omparing (2 

boundary conditions (2.15) 

number of down spins in the system changes by AA'^^ = or 1 depending on the spin of the core electron. Without 

magnetic fields the Bethe Ansatz states are highest weight in the spin SU{2), i.e. correspond to the first case. This 

results in the following expression for the exponent 



aba 



1 

2^ 



[0" - ly 



(3.16) 



From (3.15) we find that there is a discontinuity of aband{h) aX h — 1: at this point the charge bound state first 
appears leading to a jump of the exponent from 3/8 to —21/8 at u = oo (note that small negative exponents 
correspond to a 'shoulder' rather than a singularity in the absorption profileB, exponents a < — 1 will hardly lead to 
an observable feature). Large boundary potentials h — > ±cx) lead to 9^ — > '-{nc-{' \) in the strong coupling limit giving 

aband ^ -\{nl + 3r 
for finite u (see Fig. \ 

Similarly, the singularity of the abso rption intensity measured in a photoemission experiment is given by a power 
law with exponent obtained from ( 1.10 ) with AA^^ — 0: 



+ |) which is always negative. Numerical solutions of the equations show a similar behaviour 



3 1 c 2 



2e 



(3.17) 



which exhibits a jump from 5/8 to —3/8 ai h — 1 and approaches | — hnc{nc + \) aX h -^ oo for infinite coupling. 
Note that (1/2^^) varies as a function of the bulk density n^ of electrons and tJae interaction strength between 1/4 
for noninteracting fermions and 1/2 in the infinite u limit of the Hubbard modelEj, while 9'^ contains the dependence 
on the strength of the boundary potentials hi^^ (in addition to Uc and u). 

For weak boundary fields, ii-.< 1 these expressions coincide with those found in the framework of a bosonized theory 
of spin carrying electronal2l'tj provided that we identify 9'^ with the forward scattering amplitude of the core hole 
potential (see also the discussion at the beginning of Sect. y). 
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FIG. 6. X-ray edge exponents for band absorption (full line) and photoemission (dashed line) in the Hubbard model as a 
function of the boundary chemical potential h for u = \, n^ = 0.5. 

B. Absolute threshold 



Let us now consider X-ray processes which excite the system into the sector with all bound states occupied, i.e. the 
absolute threshold for absorption. Following the discussion above one has to distinguish four cases: For sufficiently 
small boundary fields (h < 1) there are no bound states, which is the situation considered in the previous section. 

For boundary fields 1 < h < u + ^Jv? -\- 1 a charge can be bound to either boundary. This changes the boundary 
phase shifts according to (|3j|) . T he computation of the finite size spectrum is complete analogeous to the case 
considere d abo ve and results in (2^). The shifts of the numbers AiV°g are now fo und t o he Q^ = ^9'^ + 1 and i 
given by (|3.11). The different boundary phase shifts modify the driving term in (3.12) to 



again 



with 



Uk) = gi°\k) + cos k Msink) 



h{x) = 2a,(,^^s){x) + i {g,^^ (|-) - G3_| (|-) } 



(3.18) 



(3.19) 



For the computation of the edge exponent from (UM we have to choose AA^^ = — 1 (the number of charges in the 
band is increased by one due to the absorption of the core electron, but at the same time two of the band electrons 
occupy the bound states in the final state). With AA^^ = as before one obtains 



(Xabs 



1 
2^ 



If 



(3.20) 



Increasing the boundary potentials such that u + y/v^ + T < h < 2u + y/ii 



1 the Bethe Ansatz state of lowest 
energy is contains both complex k and complex A leading to phase shifts ( |3.6D . As discussed above this corresponds 
to occupied charge bound states while the spinon bound states are empty. Analysing the Bethe Ansatz equations we 

find 6'^ = i^^ — 1. The fimction pc (fc) is determined by the same set of equations ( 3.12| ), (3.18) and (B.19) as above. 
The state relevant for the edge exponent is now determined by the quantum numbers AiV^ = — 1 and AA^^ — —2 
which gives again (3.20|). 
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FIG. 7. Exponents at the absolute threshold for X-ray absorption in the Hubbard model as a function of the boundary 
chemical potential h for u = 1 and several densities Uc- 



A final change in the configuration describing the absolute ground state occurs for h > 2u + y/Av^ + T {S > 2u). 
The presence of bound pairs of electrons leads to the phase shifts (3.7) in the Bethe Ans atz e quations. The quantities 
determining the edge exponents are now 0* — ^0^ , where 9^ has to be computed from (|3.12|) with 



9c{k) = g^°\k) + 2cosfc{a2s(sinfc)-a2(s_2u)(sinfc)} . 

The quantum numbers of the final state are AiV° = —3 and AiV° = —2 which gives 

1 1 

~2e 



Oiabs 



(^ + 3)' 



(3.21) 



(3.22) 



Again, the equations si mpli fy significantly in the strong c oupli ng limi t wh ere one should rescale S" by u to see the 
different regimes. Using (3.14) we can combine Eqs. (3.16), (3.2C) and (3.22) into aabs — i (l ^ ^^) where 



2 ( f h — cos TTTlj 

X = — < arctan -, 

TT V sin Trrir 



arctan 



h- Au 



(3.23) 



2 sinTrric 
Hence we find the following expression for the edge exponent of the absolute threshold in the strong coupling limit 

-i(2ric-H5)(2nc-hl) for /i< -1 
g^^,.c I ^M-"c 3) forK/K4u . (3.24) 



Oiabs 



-i(2n, + l)(2n, 
\{2n,,~2,){2n. 



7) for h:^ Au 



Since we consider the Hubbard model at less than half filling (i.e. ric < 1) this implies that a positive exponent a 
leading to a edge singularity is possible only in the intermediate regime. The corresponding numerical data for finite 
u are presented in Figure |j. 

C. Intermediate thresholds 

Finally we consider some cases where the absorption excites the system into a state in which some but not all 
bound states are occupied. First, let the final state be characterized by one antiholon and one spinon in a bound state 
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which gives rise to a singularity at an energy between the two thresholds discussed above. Such a process is possible 
for boundary potentials h > u + a/m^ + 1 (or S > u) and corresponds to a Bethe Ansatz state with a single complex 
k. A nalyzi ng th e Bet he Ansatz equations we ob tain the relation 6'^ = ^9'^. In this case 9"^ has to be computed from 



Eqs. (|3.1lD and (|3.12|) with gdk) given by ( p.lSD with 



Mx) = -a2s{x) + ^G|+i (|-) . (3.25) 



The finite siz e spe ctrum is again of the form (2.8); the quantum numbers of the relevant final state are AiV^ = 



AN^. From (|l.lOD we obtain 



«-* = i - ^ (^')' (3.26) 

for the edge exponent determining the singularity at this threshold. In the strong coupling limit wc find that a^^f 
varies between 5/8 for the empty band and —3/8 as we approach half-filling. An edge singularity can be observed for 

"c < \/i- 5 w 0.725. 

A different intermediate thershold occurs if only an antiholon is in one of the bound states. This final state is already 
possible for h > I and is parametrized by a single complex root k = ilnh for S < u and an additional complex spin 
rapidity A for S > u. Depending on h several cases have to be distinguished resulting in a edge singularity with 
exponent 



a. 



' _ 1 1 /^c^2 



o - 7T72 ((^T (3.27) 



int 2 2^2 



for S < Su (for S* > 3ii the exponent is always negative). The function fb{x) in ( 3.18 ) is now simply one half o f tha t 
in (3.19). In the strong coupling limit the edge exponent a^„j can be expressed through Uc and h using Eq. (3.15). 



In this limit a singularity in the absorption spectrum (i.e. positive exponent) can be observed for sufficiently large 
boundary potentials h > tan (^(2ric + 1)) > 1 as long as tIc < ^ but only close to /i « 4u above quarter filling. 

Note that for sufficiently strong boundary potentials the cases discussed here are only a small subset of the possible 
thresholds. Furthermore, for sufficiently strong repulsive boundary potentials, i.e. h < —1, the spectrum allows 
for holon bound states. Like in the case of the t-J model with attractive boundary chemical potentials there exist 



solutions to the Bethe equations with complex quasi momenta k = tt + i In | /i | of Eq. (3.2) that have a gap with respect 
to the absolute ground state and lead to a higher threshold in the X-ray spectrum. 

IV. CONCLUSIONS 

In this work we have determined the X-ray edge exponents in a Luttinger liquid for the case where the local 
disturbance due to the core hole leads to bound states. We used specific realizations of Luttinger liquids on the 
lattice, namely Hubbard and t-J models with integrable boundary terms. The main difference to the Fermi liquid 
case (1.1) solved in Refs. Wa is that due to spin and charge separation we find a richer structure of thresholds in the 



X-ray absorption rate associated with bound states of spinons and (anti)holons. Using Boundary Conformal Field 
Theory the exact dependence of the edge exponents on band filling and interaction strength can be extracted from 
the finite size spectra which are determined from the Bethe Ansatz solution. _. 

For weak boundary fields our results coincide with those obtained in a field theoretical treatment by Prokof'evllEl 
and Affleck and Ludwigllil if the boundary chemical potentials are fine-tuned. 

For sufficiently strong boundary fields the models considered in this paper allow for various bound states, each of 
which can lead — in principle — to a singularity in the absorption spectrum. Previous studies of these additional 
singularities,h.ave not taken into account the interaction between the particles in the bound states and those remaining 
in the bandau. This results in a simple relation between the exponents at different edges with the phase shift S{eF) at 
the Fermi surface as the only free parameter. In the systems considered here the occupation of the boundary bound 
states modifies the potential acting on the particles remaining in the bands which in turn modifies the corresponding 
phase shifts. Examining the edge exponents for the different thresholds we find that for generic values of boundary 
potentials and filling factors many of them will in fact be negative, and consequently won't lead to an observable 
singularity in the spectrum. 
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